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We present a theoretical study of exotic hybrid meson (J PC = 1 + ) production in photon-photon 
collisions where one of the photons is deeply virtual, including twist 2 and twist 3 contributions. 
We calculate the cross section of this process, which turns out to be large enough to imply sizeable 
counting rates in the present high luminosity electron-positron colliders. We emphasize the impor- 
tance of the irn decay channel for the detection of the hybrid meson candidate 7ri(1400) and calculate 
the cross section and the angular distribution for 7r n pair production in the unpolarized case. This 
angular distribution is a useful tool for disentangling the hybrid meson signal from the background. 
Finally, we calculate the single spin asymmetry associated with one initial longitudinally polarized 
lepton. 



45 
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<-| , I. INTRODUCTION 



Photon-photon collisions, with one deeply virtual photon, is an excellent tool for the study of different aspects 
of QCD. The main feature of such processes is that a QCD factorization theorem holds, which separates a hard 
partonic subprocess involving scattered photons from a distribution amplitude Q describing a transition of 
a quark-antiquark pair to a meson or a generalized distribution amplitude Q describing the transition of a 
quark-antiquark pair to two- or three-meson states ,3] . The success of this description Q with respect to recent 
LEP data |5j is indeed striking enough to propose such processes to be used as a tool for the discovery of new 
hadronic states 0. In this paper, we focus on the process where an exotic J PC = 1 ^ isotriplet hybrid meson 
H (which may be the 7Ti(1400) state [7() is created in photon-photon collisions and then decays into ir°i]. This 
exotic particle has been the subject of intensive studies for many years [8j. 

In previous papers, we have shown [9j that, contrary to naive expectations, the longitudinally polarized hybrid 
meson possesses a leading twist distribution amplitude (DA) related to the usual non-local quark-antiquark 
correlators. We were also able to estimate the normalization of this DA, which is by no means small, with 
respect to the one for usual non exotic mesons. We thus advocated that exclusive deep electroproduction 
processes will be able to copiously produce these exotic states. In this paper, we extend our analysis of (3 to 
7* 7 collisions with production of both longitudinally and transversally polarised hybrid meson. We calculate 
the hard amplitude up to the level of twist 3 and thus ignore the contributions of mass terms. Indeed, the case of 
hybrid production and its decay products, irq pair in electron-photon collisions is similar to the electron-proton 
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case, with the important distinction that no unknown generalized parton distribution enters the amplitude, 
so that the only place where non-perturbative physics enters is the final state DA or generalized distribution 
amplitudes (GDA). We emphasize the irrj pair production as a promising way for detecting the hybrid meson. 

Since we deal with the production of an isovector state, there is no mixing between quark-antiquark correlators 
and gluon-gluon correlators. Since lepton beams are easily polarized, we also consider the single spin asymmetry 
associated with the case where one of the initial lepton is longitudinally polarized while the polarizations of the 
other one are averaged over. This single spin asymmetry will turn out to give access to the phase difference of 
leading twist and twist 3 components of the final state GDA. 

Finally let us note that in the case of the electron-photon collisions producing a ttt] pair, i.e. a state with the 
positive charge parity C = +, the Bremsstrahlung contribution does not exist. 



II. VACUUM-TO-HYBRID MESON, VACUUM-TO-HADRONS MATRIX ELEMENTS AND 

THEIR PROPERTIES 



The hybrid meson Distribution Amplitude 

Let us first consider relevant vacuum-to-hybrid meson matrix elements which enter in exclusive hard am- 
plitudes. We suppose that the hybrid meson is in the isotriplet state with J PC = 1 h quantum numbers. A 
meson with such quantum numbers should be built of quarks and gluon, that is one should work beyond the 
quark-antiquark model. We have shown || that in the case of the longitudinal hybrid meson polarization the 
leading contribution to a hard amplitude comes from the non-local quark-antiquark correlators. The non-local 
character of such correlators leads to the gluonic degrees of freedom included in the gauge-invariant link. 
Let us form the light-cone basis adapted for 7*7 — > H and 7*7 — > irij processes. We introduce the "large" and 
"small" vectors as 

n* = (A,O r ,A), n = (-Lo r ,~), n*-n = l, (1) 

respectively, and express the photons and hybrid meson momenta via the Sudakov decomposition: 

. Q 2 , Ml + Q 2 M% 
q = n - —n, q = n, p = n + —f-n. (2) 

We define the transverse tensor g^ v — g^ v — n*n v — n^n*. 

Using the results of seminal studies on the DAs of vector mesons ^(J , we write (u = 1 — u): 



(3) 



(H(p,X)\^(z)^[z;-z]^(-z)\0} = 
1 

f H M H Pll e w -n J due iiu -^ p - z <t>?(u)+e { *l \ due i{v -*> z (t>% \u) 


for the vector correlator, and 

(H(p,X)\i,(z) lll75 [z;-z}i>(-z)\0) = if H M H E n ^ )pn \ du e *(«-«)f» ^(u) M ) 







for the axial correlator. We use the following short notation : £ s kmt — e MiM2A'3M4' s Mi^M2 m A'3^4- ^ n ecms © an d 
|Hl. the polarization vector describes the spin state of the hybrid meson. Due to C-charge invariance, the 
symmetry properties of DAs are manifested in the following relations: 

4>f (u) = -4>f (1 - u), 0?(u) = -tff(l - u), ^(«) = ^(l-u). (5) 

Compared to the p meson matrix elements [lo|. one can see that the corresponding DAs for the exotic hybrid 
meson has different symmetry properties. 
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The leading twist longitudinally polarized hybrid meson DA has been discussed in |9| and shown to be 
asymptotically equal to 



<${u) =30u(l-u)(l-2tt) 



(G) 



and fjj to be of the order of 50 MeV. In the present work, we will not include any QCD evolution effects in 
GDA. 

In order to be able to estimate twist 3 contributions, we have to build a model for the twist 3 D As 6 ^ and 
4>\. We will not innovate on this point but restrict to the Wandzura-Wilczek parts which read [l0| : 



yielding for our choice of <j>^ 



dv 



v-1 



v — 1 



(7) 



±WW r \ 
03 H = 



--(l~2uf 
|[l + 6«(«-l)] 



(8) 



The 7r?7 Generalized Distribution Amplitude 

We now come to a discussion of the irn GDA which enters the amplitude, if we try to detect the hybrid meson 
through its nrj decay mode, which may be the dominant one for the 7Ti(1400) state. Using the results of Ref. 
Hal we have 



(7T°( P7r ) V (p v )\4>(-z)^[~z; z]iP(z)\0) = 

l l 

Pnn J due^-^^ („, C, ml n ) + A^ J due^^^ («, C, m^) 

o o 

for the vector correlator, and 

l 

o 

for the axial correlator. In @ and 1|10|) the total momentum and relative momentum of irr/ pair are 



in 
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The ttt] leading twist GDA has been discussed in [jj in the J PC = 1 + channel as 



^P\u, C, mi n ) = 30u(l - u)(l - 2u) B 11 (m' 7rri )Pi(c 



(9) 



(10) 



2 2 
Tfl — 777 

P* v = Ptt +P V =n* + -ZL n , A nn =p w - Pri = (2t-l+ ^ 2 V )P WV + (1 - 2()ml v n + A^. (11) 



(12) 



with cosQ = (2( — l)//3, (3 = A(?n^. , m^,m^)/m^ , and the coefficient function Bn(m% ) related to the Breit- 
Wigner amplitude when m^L is in the vicinity of Mfj as 



-Bii(m^) 



5 gHjnjHMH^ 

3 Ml - ml v - iY H M H ' 



(13) 
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The formalism of GDA includes a background in a natural way, i.e. the contribution of final states uncorrelated 
with the main signal of hybrid meson exchange. We will model this background as a J = contribution (i.e. 
(■—independent) without structure in m^, and with the asymptotic u— dependence (i.e. u(l — u)(2u — 1)). We 
know nothing of its phase but that it should be quite constant in the considered region. For simplicity, and 
contrary to the analysis of Ref. ||, we do not include any 3 = 2 component to the GDA. Thus our model for 
the leading twist irrj GDA reads 



$i'""(u»C, n &7 ) =30u(l-u)(l-2u)(Ke ia + B u (ml Jcos 



(14) 



Note that the GDAs $3 (u, £, Tri^ n ), 9^ (u, £, m^.„) are twist 3 functions; they have a part which comes from 
the Wandzura-Wilczek (WW or kinematical) twist 3 and another part which is usually called genuinely twist 3 
(see, for instance [l2j). The WW part reads 



dy- 



dy- 



(15) 
(16) 



which yields with our choice of the leading twist terms: 



$3^( M ,C)=5(l-2 U ) 3 ^ig^ 

*y («, C) = -5 ((1 - 2uf - 2u(l - „)) E}^A 



(17) 



III. AMPLITUDES OF 7*7 -> H AND 7*7 -> ivq PROCESSES 



The gauge invariant expression for the ~/{q')"f*(q) — * H(p) amplitude reads 



1 flU 



e 2 {Ql-Ql) 
V2 



1 T (X) \ (P + 9% a 

-g^e w ■ nAi + — 2 ^A 2 



with Q u = 2/3 and Qd = —1/3, and where the following short notations have been introduced: 



(18) 



1 1 

Ai = j ' duE-(u)$i(u), A 2 = j du^E-(u)$ 3 (u) - E + (u)<S> A (u) 


$1,3,4") = fnM HC f>f 3 A (u), E ± {u) = — !— ± -. 

1 — u u 



(19) 



Such an amplitude allows calculating the production cross section and the hclicity density matrix of the hybrid 
meson. The knowledge of this helicity matrix leads to a definite angular distribution for any particular decay 
channel. In particular, for the process leading to the two meson ■nrj final state the straightforward generalization 

:eads 

e 2 (Ql-Ql) 



of the amplitude derived in 



ie_pr 

.12] 



r Fl1 — * 7T7 1 



V2 



-a 1 A y 



(7TT7J 



Q 2 



(20) 
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where, following to the notations l|19|l . we introduce: 



= J duE^(u)^ v \u) = -10 [Ke la +B 11 (ml n )cose] 



l 

4^ - J du^(u)^\u)~E + (u)^r\^) 



5 

"3/3 



(21) 



Note that the amplitudes l|18(l and l|20|) satisfy the gauge invariance condition: q^T^ = q'JT^ = provided we 
neglect terms proportional to the square of the meson mass, which is quite natural in the light-cone formalism 
in the twist 3 approximation. 



IV. CROSS SECTIONS 



The kinematics of 7*7 — > nrj process is illustrated in Fig. ^ We can now calculate the cross sections for the 

, y 




Figure 1: Kinematics of the process e 7 — > e n r\ . 

process 7*7 — ► Hybrid meson. As for the usual treatment of a pseudoscalar meson it may be expressed in 
terms of the transition form factor Fn-y, which scales like 1/Q 2 up to logarithmic corrections due to the QCD 
evolution of the DA (which we consistently ignore in this work). For an easy comparison with the well measured 
cases, we define R as the ratio of squared amplitudes for unpolarized photons, 



R 



r^(77* -> Bjr*^rf -> h) 



(22) 



As shown on Fig. by the solid line, the twist 2 transition form factor is sizeable in the hybrid case, of the 
order of the corresponding quantity for 7r° or r\ meson production. We are thus confident that a good detector 
as the ones existing in the present e + e~ colliders will be able to detect the hybrid signal if one of the decay 
channel has a fairly large branching ratio. To estimate the twist 3 effects, we now approximate the twist 3 part 
of the hybrid distribution amplitude in the Wandzura Wilczek way [Tlj as explained in Section II. The resulting 
ratio R with the twist 2 and twist 3 contributions to the hybrid transition form factor (but only the twist 2 
contributions to the tt° case) is shown on Fig. [3 by the dashed line. The twist 3 contribution is negative and 
of the order of 20% when Q ss 1 GeV, but quite negligible when Q > 3 GeV. Let us stress that this estimate 
is in the WW approximation which is by no means a proven result. Calculating the order of magnitude of the 
genuine twist 3 hybrid DA is very model dependent, and to our knowledge, no estimate exists in the literature. 

The cross section for the e + e~ process 



e(fci) + e(h) -> e(fe) + e{h) + H(p). 



(23) 



for unpolarized lepton beams is easily obtained when including the leptonic parts. Note that the positive C 
parity of the hybrid meson does not allow any contribution from a Bremstrahlung process. Specifying a positive 
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Figure 2: The ratio R(Q 2 ) of the squared amplitudes for H and 7r° production in 7*7 collisions at leading twist and 
zeroth order in a s (solid line) and including twist 3 contributions in the numerator (dashed line). 



C parity two body decay channel, like 7r°77, we get the differential cross section for the complete process (see 
Fig.EJ 



i{ki) + e(h) -> e(k 2 ) + e(/ 2 ) + ^{p-k) + viPv)^ 



(24) 



where we average (sum) over the initial (final) lepton polarizations and use the equivalent photon approximation 
(EPA) for the leptonic part connected to the almost real photon, 



da,, 



a 1 1 



(l-z 2 ) 



dQ 2 dW 2 d cos 6 dip dx2 it X2 



■In 



Qmax{ X ^) 



d(Jr 



dQ 2 dW 2 d cos 6 d^ 1 



(25) 



where W = m„, and Q' m i n an d Q'max are the minimal and maximal virtuality of the photon q' , respectively. 
For a given ee collider energy, the variables X2 = q'p/hp = s e -y/s ee and y = qq'/kiq' are not independent at 
fixed Q 2 and W 2 , since yx2 — {Q 2 + W 2 )/s ee . tp is defined as being the angle between leptonic and hadronic 
planes. The lower kinematical limit Q' min = x\m 2 /(l — X2) is determined by the electron mass m e , whereas 
Q'max depends on experimental cuts. Without a precise knowledge of these cuts, we will present results for the 
ej — » e ttt] process. Its cross section reads 

a 3 f3 1 I 
2j//? 2 W 2 



(Z(7 e .-y— > e tj-^ 

dQ 2 dW 2 d cosO dip 

/i + (i-y) 2 



(26) 
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(71-77) 1 2 



■sin 2 0|4 



Qy 2 



cos(^sin0Re(4"' 7) 4 X " ) *) 



V 4y 2 1 1 1 Q 2 y 2 

To analyze the cross section, let us first define the integrated over angles 9 and <p cross section 
da(Q 2 ,W 2 )/dQ 2 dW 2 



da(Q 2 , W 2 



d cos 9 dtp- 



da,, 



25a 3 em {3 l + (l-y) 2 



72s 2 Q 2 



dQ 2 dW 2 J r dQ 2 dW 2 d cos 9 dip 

where we restrict ourselves to the twist 2 contribution and 

5 



K 2 



1 



a 2 {W) b 2 {W) 



a{W) = -g H7Iv f H PM H , 6(W0 



1 



(27) 



(28) 



V(W a -Jlf£)2 + I*M£ 

In Eq. I|28|) the coupling constant gu-KT) is related to the branching ratio BR = Th^^/Yh of the hybrid meson 

by 

48tt 



2 



(3 3 M H 



r H br 



(29) 
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Since nothing is known about the value of the H — > nr/ BR, we choose in our numerical analysis as a reasonable 
estimate BR = 20%. Let us also note that simultaneous rescaling of the background magnitude K and the BR 
by a coefficient c amounts to a rescaling of the cross sections by c 2 , leaving the angular distribution W(9,<p) 
(Eq. (jSOJ) unchanged. 
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Figure 3: The differential cross-section for ivr] pair production as a function of W for Q — 3 GeV, y = 0.3, for different 
background magnitudes K = (solid curve), 0.5 (dashed curve) and 1 (dotted curve) 




Figure 4: The differential cross-section for nr/ pair production as a function of W for Q = 5 GeV, y = 0.3, for different 
background magnitudes K = (solid curve), 0.5 (dashed curve) and 1 (dotted curve) 

This cross section does not depend on the phase a appearing in our model of the 71-77 — GDA, H4j |. The 
plots on Figs. |3] and 0] show its W dependence for different magnitudes K of the assumed background, for 
respectively Q — 3 GeV and Q — 5 GeV. We observe that the presence of the hybrid peak around W = 1.4 
GeV is hardly affected when changing the magnitude K of the background. The magnitude of this signal is 
comparable with what has been achieved by the L3 experiment at LEP 0. As expected, the comparison of 
Fig. Eland Fig. 0] shows that the magnitude of the signal decreases when Q increases. However, this decrease is 
not so dramatic due to the scaling behaviour of the amplitude from twist two dominance. We thus expect that 
the Q dependence may be experimentally studied up to a few GeV. 
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A more detailed test of the nature of an eventual signal may be accessed by a study of the angular distribution 
of the irr] final state. Using l(2"7jl and (|2"T|) supplemented by l(r71l2"tj|) and l|2"Tl) . we obtain 

fda(Q 2 ,W 2 )Y da e 7 ^ e nv 1 . 2 



\ dQ 2 dW 2 J dQ 2 dW 2 d cos 9 dip \-k 

(30) 

with the functions A(W), B(W), C(W), D(W,Q) and E(W,Q): 



A 



K 2 



K 2 



Ln2 



b 2 



B = 



2Kab cos (7 — a) 



K 2 + \a 2 b 2 



C = 



a 2 b 2 



K 2 



D = 



C 



Q 3[l + (l-y) 2 ] 
and where 7 is the phase of the Breit-Wigner form, i.c 



Q 3[1 + (1 - y) 2 ] 



cos 7 = (M| - W 2 ) b(W), sin 7 = b(W) T H M H 



(31) 



(32) 



In (|31() we restrict ourselves to the twist 2 contributions for the coefficients A(W), B(W), C(W), 
which is why A, B, C are Q independent. Due to the normalization of the angular distribution W(0,ip), 
J d cos 9 d ipW(9, <p) = 1, the functions A and C are not independent but they satisfy the relation A + C/3 = 1. 
The function A is displayed in Fig. Its shape is dictated by the inverse shape of the integrated cross-section 
shown in Fig. |3| as seen from (|31|) and (|27|) . This is intimately related to our main assumption that the 
background is described by a J = contribution. 
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Figure 5: The A component of the angular distribution function as a function of the nr/ mass W for Q — 3 GeV, y = 0.3, 
for different background magnitudes K = 0.1 (solid curve), 0.5 (dashed curve) and 1 (dotted curve) 

The B coefficient is displayed in Figs. and [7] for a background magnitude K = 1 and y = 0.3. It measures 
the interference between the background and the hybrid signal. It is thus quite dependent on the value of the 
phase of the background, but its W dependence always reveals a dramatic change around the mass of the hybrid 
meson. One may use it to determine more precisely this mass. 

Since D and E vanish at the twist 2 level, this part of the angular distribution of the final mesons is sensitive 
to the strength of the twist 3 amplitude. We calculate them in the Wandzura-Wilczek approximation described 
in Section II. The coefficient D is independent of the phase of the background. We show on Fig. [Slits behaviour 
when the magnitude of the background varies, for Q = 3 GeV and y = 0.3. On Fig. [5] we show its W dependence 
for K = 1 and y = .3 and for three different values of Q. Because of its proportionality with C, and thus its 
relation to A, it is strongly peaked around the hybrid mass. It will only be measurable at fairly small values of 
Q. E(W, Q) depends on the background phase a, as does B, and on the magnitude K. We show in Figs. ITU1 and 
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Figure 6: The B component of the angular distribution function as a function of the irr] mass W and of the background 
phase a for the background magnitude K = 1, with Q = 3 GeV and y — 0.3. 

B 




W (GeV) 

Figure 7: The B component of the angular distribution function as a function of the irr) mass W for Q = 3 GeV, y = 0.3 
and K = 1, for different background phases a = (solid curve), 7r/4 (dashed curve) and n/2 (dotted curve). 

Illlits W dependence for a = and K = 0.5 and K = 1 respectively. This behaviour is similar to the one of B 
displayed in Fig. [7| but its magnitude is much smaller when Q is greater than 5 GeV. 



V. SINGLE SPIN ASYMMETRY 



We consider now the exclusive process where a longitudinally polarized lepton (with helicity h) scatters on 
an unpolarized photon to produce the lepton and the hybrid meson detected through its decay into a ixr\ pair. 
Such a process allows to define an asymmetry which is zero at the leading twist level but receives contributions 
from the interference of twist 2 and twist 3 amplitudes. This asymmetry is related to the azimuthal angular 
dependence of the polarized cross section and it is defined as 
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Figure 10: The E component of the angular distribution function as a function of the 7r?7 mass W for y = 0.3 and 
K — 0.5, for different values of Q =3 GeV (solid curve), 5 GeV (dashed curve) and 7 GeV (dotted curve). 
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Figure 11: The E component of the angular distribution function as a function of the nr/ mass W for y = 0.3 and K = 1, 
for different values of Q =3 GeV (solid curve), 5 GeV (dashed curve) and 7 GeV (dotted curve). 



-,2 W 2. \ _ / dcos9 cm (da^ - da^) 



Ms^w-M = J IdcoseUd ^ )+da{ ^ , (33) 

where we denote by dfj^~*^ the differential cross section daijlj^ri /dW 2 dQ 2 d cos 9 cm dip. 

Note that the denominator (we restrict ourselves to the dominant twist 2 component) can be expressed 
through the unpolarized differential cross section defined in Eq. I|27|l 



2tt 



dtp J dco S cms (^H^H^a ^y , (34 ) 

-1 

The asymmetry (|33[1 reads: 

7T 

/ e$ cm sin0 cm .A/ r (0 cm , Q, W, ip) 

Ai (s e7 ,Q 2 ,W 2 ;ip) = ^ (35) 

2 / d^cm sin 9 cm D(9 cm ) 



where 

A^ cm , 0, W, <p) = ^e kiqlk ^ MAf'^W), p ( , cm ) = + W)2 |A^| 2 . (36) 



The crucial dynamical quantity probed by this asymmetry is thus Im(A^ 7rl, ' ) A* 2 ), It depends on the phase 
structure of the amplitudes, and may be written as 

t fA ("v) A *{*r,h 50 a(W) b(W) K . 

Im(Ai "A^ ") = - sm(a-7j , (37) 

where a(W), b(W) and 7 are defined by equations (|28|l and l|32l) . This quantity depends much on the unknown 
background phase a. On Fig. 1121 , we present our result with the choice a = 0. The resulting asymmetry is 
sizeable and should be measurable. 

In order to get more statistics, one may define an integrated asymmetry : 

2tt 

2-7T J d(cos6 cm ) J d(p sin^der^ — da^~>) 

A 2 (s ei ,Q 2 ,W 2 ) = °- . (38) 

J d(cos8 cm ) J dip (da(~*) + da(*~)) 
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Figure 12: The single spin asymmetry Ai as function of ip = (0, 2-k). Values of parameters: W = 1.4 GeV, Q 2 = 
5.0 GeV 2 , s e7 = 10 GeV 2 , a = 0. The solid line corresponds to K = 0.8, the short-dashed line to K — 1.0, the long- 
dashed line to K = 0.5. 

which may be used as a probe of the wrj invariant mass dependence. We show on Figs. 1131 and 1141 the resulting 
W dependence for two choices of the background magnitude K. 
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Figure 13: The integrated single spin asymmetry Ai as a function of the n-q invariant mass. Values of parameters: 
Q 2 = 5.0GeV 2 ,s e7 = 10 GeV 2 , K — 1.0 , the solid line corresponds to a = 0, the short-dashed line to a = 7r/4, the 
long-dashed line to a = 7r/2. 



CONCLUSION 

This theoretical study shows that if a hybrid meson exists with J PC = 1 ^ around 1.4 GeV and with a 
sizeable branching ratio to ir — i], much can be learned about it from the experimental observation of 7*7 
reactions and the precise study of the tt — r\ final state. The magnitude of the cross section that we obtain in 
our model of the 7r — 77 Generalized Distribution Amplitude indicates that present detectors at current e + e~ 
colliders are able to get good statistics on these reactions, provided the tagging procedure is efficient. 
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Figure 14: The integrated single spin asymmetry Ai as a function of the n-q invariant mass. Values of parameters: 
Q 2 = 5.0GeV 2 ,s e7 = lOGeV 2 , K = 0.5 , the solid line corresponds to a — 0, the short-dashed line to a — it/A, the 
long-dashed line to a — tt/2. 
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